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We discuss the current status of our calculation of the physics of % and K mesons using three 
dynamical flavors of improved staggered quarks. This year, we have a new ensemble with a lattice 
spacing of 0.06 fm and a light sea mass of 0.2m s , as well as significant increases in statistics at 
several coarser lattice spacings and/or heavier sea masses. Results for decay constants, quark 
masses, low energy constants, condensates, and V IIS are presented. 
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Table 1: Lattice parameters. The lattice spacings are the "nominal" scales (see text). The n and p referred 
to are those formed out of the sea quarks for each lattice; valence quark masses however go down to the 
lightest sea-quark values in the table. 



We are using improved staggered quarks [[T]] with Nf = 3 dynamical flavors (both unquenched 
("full") QCD and partially quenched) to study the physics of light pseudoscalars (71, K). Since 
our original published work [g], we have continued to add data sets with lighter sea quark masses 
and/or finer lattice spacings and to improve the analysis. This is the latest in a series of periodic 
updates [Q, H]. We concentrate here on those aspects that have changed since last year. 

Table g gives the parameters of our lattices. The quantities m' s and m' = m' u = m' d denote the 
values of sea quark masses chosen in each run. (The corresponding masses without the primes, 
e.g., m s and m = (m u + nid)/2, are the physical values.) 

The aw 0.06 fm lattice with masses 0.0036 / 0.018 is a new ensemble this year, as is (for this 
analysis) the large- volume a^0. 12 fm lattice with masses 0.01 / 0.05 and spatial size 28 3 . The 
numbers of configurations for the am 0.06 fm lattice with masses 0.0072 / 0.018 and for several of 
the aw 0.1 2 fm lattices have almost doubled since last year. Running on an aw 0.06 fm lattice with 
m' = 0.1m', (masses 0.0018 / 0.018) has recently begun but is not included here. 

On each ensemble, we determine r\/a, where r\ (ml \m' s ,g 2 ) []|] is a length scale from the static 
quark potential, similar to ro [0]. The quantity rf 1 ^, defined as the continuum r\ at physical quark 
masses (m, m s ) may be determined from the r\ja values and the T 2S-1S splitting We obtain 
rf ys = 0.318(7) fin flUl. 

For generic chiral and continuum extrapolations, it is convenient to define the lattice scale by 
a = r\ y& / (r\(m! ,m' s ,g 2 ) /a). We call this the "nominal" scale-setting procedure. In choosing the 
input lattice coupling g 2 , we kept r\ja fixed as m! and m' s changed over a given set of ensembles 
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{e.g., the a w 0.12 fm ensembles). Thus, up to tuning errors, each ensemble grouped within a box 
in Table [I] has the same nominal scale. However, fixing the scale this way is not completely correct 
for applying chiral perturbation theory (ZPT) to quantities such as f n since r\ has some (small, but 
physical) dependence on the dynamical quark masses that is not included in ZPT 

A mass-independent procedure to set the scale is preferable [Q]. A convenient procedure is to 
replace ri(m',m' s ,g 2 )/a by ri(m,m s ,g 2 )/a, where the value of r\(m,m s ,g 2 )/a at physical masses 
m,m s is obtained by a smooth interpolation/extrapolation from ri(m' \m' s ,g 2 )/a. We tried this mass- 
independent scheme in Ref. but the differences with the nominal approach were smaller than 
other systematic errors for all quantities. With better data, we now find significant differences in a 
few low energy constants (LECs). In addition, the mass-independent scheme tends to have better 
confidence levels in our XPT fits. Therefore we use this scheme exclusively here. 

As in Refs. [Q, ^, |[] we fit the partially quenched (PQ) lattice data to rooted staggered chiral 
perturbation theory (rSZPT) forms [10, [H], 12]. We always fit multiple lattice spacings, and both 
masses and decay constants, simultaneously. To determine the LO and NLO LECs and chiral- 
limit quantities, we fit to the low quark-mass region, and omit the aw 0.15 fm lattices, where taste 
violations are large. Denoting the valence quark masses in the mesons by m x and m y , the low- 
mass cuts are: am x + am y < 0.39 am s (at aw 0.12 fm); am x + am y <0.51 am s (at aw 0.09 fm); and 
am x + am y < 0.56 am s (at aw 0.06 fm). We can tolerate a higher cutoff at smaller lattice spacing 
because the taste violations, and hence the masses of non-Goldstone pions, are smaller. In these fits, 
we also cut on sea-quark mass and remove the aw0.12 fm sets with masses 0.03/0.05, 0.02/0.05, 
and 0.03/0.03. Because the statistical errors are so small, we still need to add in the NNLO analytic 
terms to the complete NLO forms in order to get good fits 

For interpolation around m s , we must include higher quark masses. Once LO and NLO pa- 
rameters are determined, we fix them (up to statistical errors) and fit to all sea mass sets, all lattice 
spacings, and valence masses m x + m y < 1.2 m s . We now also need to add in NNNLO analytic 
terms to get good fits. These NNNLO fits are used for central values of f n , f% and quark masses. 

Figure [I] shows results for the squared pseudoscalar masses as a function of quark mass. "Pi- 
ons" have valence masses m y = m x ; while "kaons" have m y held fixed at various (arbitrary) values 
while m x varies. The fit is to the full quark-mass range and uses NNNLO terms. 

For the pions, the relative values of the results on various lattices is determined largely by the 
relation between the simulation strange mass m' s in the sea and the physical mass m s . For example, 
m'Jm s is largest for the aw 0.12 fm lattices, which makes the slope of the pion data greatest for 
these lattices. For kaons, the biggest effect is simply the choice of the values of the fixed valence 
mass rriy, typically chosen to be various fixed fractions of m' s . 

Extrapolating to the continuum and setting valence and sea quark masses equal, we get the 
dashed red lines; m' s has been adjusted so that both the kaon and the pion hit their physical values 
at the same value of m x . This gives the physical quark masses m and m s (after renormalization). 

Note that the fit lines in Fig. are remarkably straight on this scale. To see curvature coming 
from the NLO chiral logs as well as the analytic higher order terms, we plot m 2 /(m x + m y ) in Fig. || 
(left). As the lattice spacing decreases, the PQ log at small mass becomes more evident. At larger 
lattice spacing, the PQ log is largely washed out by staggered taste violations. The continuum 
dashed red line has sea and valence masses equal, so no PQ log is expected. 
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Figure 1: Comparison of NNNLO fit to partially-quenched squared meson masses. For clarity only the 
lightest sea-quark ensemble for each lattice spacing is shown. 
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Figure 2: Data for m^/(m x + m y ) is plotted at left; while that for f K is plotted at right. This is the same fit 
as in Fig. 1 . All sea-quark ensembles are represented, but only "pion" points (m x = m y ) are shown. 



In Fig. |2| (right), we show the behavior of the decay constant. Extrapolating to the continuum, 
setting m' s = m s , and setting light valence and sea masses equal gives the dashed red line. The final 
result for f„ after extrapolating m x ,m y —>-rh is marked by a +. The experimental result is indicated 
by the o; it comes from the 7i + -> /i+ + decay width and V ud = 0.97377(27) Q. 

All points and fit lines above have been corrected for finite volume effects using the rSZPT 
forms at one loop. However, it is known [14] that finite volume effects coming from higher orders 
in XPT can be a large (~ 50%) correction to the one-loop effects in the current ranges of quark 
mass and volume. We therefore study this issue directly by comparing results on the spatial size 
20 3 and 28 3 a«0.12 fm ensembles with am' = 0.01, am' s = 0.05. These lattices have spatial length 
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1-loop % diff. 




1.4(2)% 


1.6(2)% 


1.1% 
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0.4(3) % 


0.4(3) % 


0.3% 


{am n ) 2 


-1.0(4)% 


-1.2(4)% 


-0.9% 


(am K ) 2 


-0.4(2) % 


-0.4(2) % 


-0.2% 



Table 2: Finite volume effects. The second column is the % difference between the values on the 28 3 and the 
20 3 lattice. In the third column, we "boost" the difference, to take into account the (small) further difference 
between 28 3 and infinite volume. (We use one-loop results to make the adjustment.) The last column shows 
the % difference between 20 3 and infinite volume as predicted by one-loop rSZPT. "Pion" quantities are 
from lattice valence masses 0.005, 0.005; "kaon" quantities are from lattice valence masses 0.005, 0.04. 

2.4 fm and 3.4 fm, respectively. The comparison is shown in Table ||. 



As expected from Ref. Q14[], the true finite volume effects are larger than those predicted at 
one-loop, although only for f n are the relative errors small enough to make the comparison unam- 
biguous. We define the "residual finite volume effect" on the 20 3 lattice as that effect not taken into 
account at one-loop, i.e., the difference between columns three and four in Table ||[ In practice, the 



20 3 , 



am 



0.01, am' = 0.05, a « 0.12 fm lattice is close to the worst case in our data set, since the 



volumes both at the lightest sea quark masses and at the finest lattice spacings are larger. 

Judging by the one-loop results, we expect the overall (and hence residual) finite volume ef- 
fects closest to the chiral and continuum limits in our data set to be about half those seen in the 
table. We therefore correct our data by 1/2 the residual finite volume effects from Table ||, and take 
the full size of the correction as a systematic error. We note that the size of the error determined 
here is very similar to that estimated by us previously [§] using Ref. [14]. 

The fact that we can get good fits to the forms predicted by rSZPT (and not to those of con- 
tinuum XVT [§]) is an overall test of staggered chiral perturbation theory, including the "replica 
trick" to represent rooting. As a more focused test of the replica trick in rSZPT, we allow n r , the 
number of replicas per staggered flavor, to be a free fit parameter. If rSZPT is correct, we should 
find n r = 1/4. On the low-mass data set described above, we obtain n r = 0.28(2)(3), where the 
errors are statistical and systematic (describing the variation over details of the chiral fits). While 
the ability of rSZPT to describe rooted staggered data cannot prove the correctness of the rooting 
trick itself, it does indicate that no problems occur in the chiral sector of the rooted theory [|I2|]. 
This is because rSZPT reproduces continuum ZPT in the limit a — > 0. 

Using r\ = 0.318(7) fm from T splittings, we obtain (still preliminary) 

f x = 128.3 ±0.5 IfsMeV 
f K = 154.3 ± 0.4 IjjMeV 
fit/ fx = 1.202(3X^1) , 

where the errors are from statistics and lattice systematics. These results are consistent with our 
previous results [Q], with 20-30% smaller errrors. Our value for f n is consistent with the experi- 
mental result, /£ xpt = 130.7 ±0.1 ± 0.36 MeV 

Instead of setting the scale from T splittings, we can set the scale from f K itself, which gives 
smaller errors for %-K quantities. Note that even dimensionless quantities can change with the new 
scale, due to changes in physical quark masses. We then obtain (preliminary): 
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h = 156.5 ±0.4 t\n MeV f K /f n = 1. 197(3) (t 13 6 ) 

hlh = 1.052(2)(±f) (uu) 2 = -(278(l)(jf )(5) MeV) 3 



f n lh = 1.21(5)(+ 13 ) (uu) 3 = -(242(9)(+ *)(4) MeV) 

-38 \ 
-15 J 



/ 2 // 3 = 1.15(5)C 13 ) (uu) 2 /{uu) 3 = 1.52(17)(+f ^ 



2U-U = 0.4(1) (jf) 2L 8 -L 5 = -0.1(1)(1) 

L 4 = 0.4(3)(t 3 ) L 5 =2.2(2)(±?) 

L 6 = 0.4(2)(^) L 8 = 1.0(l)(l) 

m s = 88(0)(3)(4)(0) MeV m = 3.2(0)(1)(2)(0) MeV 

m u = 1.9(0)(1)(1)(1) MeV m d = 4.6(0)(2)(2)(1) MeV 

m s /m = 27.2(1) (3) (0)(0) m u /m d = 0.42(0) (1)(0) (4) . 

The errors are statistical, lattice-systematic, perturbative (for masses and condensates; from two- 
loop perturbation theory [|15|]) and electromagnetic (for masses; from continuum estimates). f 2 (f$) 
represents the three-flavor decay constant in the two (three) flavor chiral limit, and (uu}2 {{uu)^) is 
the corresponding condensate. The low energy constants L ; are in units of 10~ 3 and are evaluated 
at chiral scale m^, the condensates and masses are in the MS scheme at scale 2GeV. 
We also obtain 

ri =0.3108(15) (^) fin, 

which is 1-c lower (and with somewhat smaller errors) than the value from the T system. There 
is a 2-a conflict between our r\ result from f n and the HPQCD Collaboration [^] value from T 
splittings, r\ = 0.32 1 (5) fm. If instead we compare to our own evaluation of r\ from the T spectrum, 
T\ = 0.318(7) [Q], the difference is only 1-a. We emphasize, however, that the evaluations of 
r\ from the T splittings both by us and by the HPQCD Collaboration use the same lattice data: 
HPQCD T splittings [0] and MILC values of r\/a The difference is only in how we extrapolate 
to the physical point and estimate the systematic error. Our result is consistent, though, with the 
(Nf = 2) result from the ETM Collaboration ro = 0.454(7) fm. Converting from tq to r\ using 
the ratio ro/n = 1.46(1)(2) (from Ref. Q|, adjusted for the slight difference between Nf = 3 and 
Nf = 2), this gives r x = 0.311(7) fm. 



Together with the experimental result for the kaon leptonic branching fraction [17], our result 



for fK/fn implies \V US \ = 0.2246(^^3), which is consistent with (and competitive with) the world- 
average value \V US \ = 0.2257(21) [ |l3| ] coming from semileptonic ^-decay coupled with non-lattice 
theory. 

The change in the perturbative mass renormalization constant Z m from one to two loops ac- 
counts for almost all of the difference between the mass values quoted here and those in Ref. Jl8|,|2|]. 
A non-perturbative evaluation of Z,„ is in progress. 

We stress that our extraction of the L, uses fits that include (analytic) NNLO terms. Therefore, 
a comparison to other evaluations, either phenomenological or on the lattice, that stop at NLO 
terms is problematic. Indeed, NNLO terms of "natural size" in ZPT can produce changes in the L, 
(relative to a pure NLO evaluation) that are as large as, or even somewhat larger than, our current 
systematic errors. This is confirmed by NLO fits to our data. Such fits have very poor confidence 
levels, however, which is why we do not include them in the analysis. 
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The SU(2) L x SU(2) R LECs l 3 , l 4 that are extracted from our SU(3) L x SU(3) R results 
using one-loop (NLO) formulae are therefore quite sensitive to the NNLO terms, particularly for 
T3. An NLO fit, on the other hand, gives 1$ = 2.85(7) (statistical errors only), which is comparable 
to the results from groups [ 16, 20] performing two-flavor simulations with NLO SU (2)l x SU (2)r 
fits. Indeed, this must be true, because the m\ data are so linear (see Fig. []]), which requires Z3 to 
have roughly this value [jig]. Alternative fits using SU (2)i x SU (2)« rSZPT are in progress. Since 
the strange sea-quark is omitted from the chiral theory, the approach should make possible good 
NLO fits on light-mass data, and thereby bypass this issue. Inclusion of two-loop (continuum) 
chiral logs [21] in SU (3)l x SU (3)# fits is also in progress. 

This work is supported in part by the US DoE and NSE Computations were performed at the 
NSF Teragrid, NERSC, and USQCD centers, and at computer centers at the University of Arizona, 
the University of California at Santa Barbara, Indiana University, and the University of Utah. 
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